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Abstract 

The Ising model with ferromagnetic interactions that decay as is ana- 
lyzed in the non-extensive regime < a < d, where the thermodynamic limit 
is not defined. In order to study the asymptotic properties of the model in 
the — > oo limit {N being the number of spins) we propose a generalization 
of the Curie- Weiss model, for which the N ^ oo limit is well defined for 
all a > 0. We conjecture that mean field theory is exact in the last model 
for all < a < d. This conjecture is supported by Monte Carlo heat bath 
simulations in the d = 1 case. Moreover, we confirm a recently conjectured 
scaling (TsalhsEl) which allows for a unification of extensive (a > d) and 

non-extensive {0 < a < d) regimes. 
75.10.Hk, 05.50.+q., 64.60.Fr 
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It has been known for a long time that systems with long-range microscopic interactions 
can exhibit non-extensive behaviour (see Refs.[l-3], among many others, and references 
therein). In other words, if the effective range of the interactions between the constituent 
particles decays slowly enough with the distance, the free energy F = —(3 In Z, with Z = 
Tr exp{—(3 H) {H being the Hamiltonian of the system and (3 = ksT) will grow /aster than 
the number N of microscopic elements when N oo and the so-called thermodynamic limit 
will be not defined. 

Besides their fundamental theoretical interest in physics, microscopic models with long- 
range interactions which decay slowly are of interest nowadays, in view of their relationship 
with neural systems modeling^, where far away localized neurons interact through an action 
potential that decays slowly along the axon. Another related problems, are spin sytems 
with RKKY like interactions, which are present in spin glasses!, critical phenomena in 
highly ionic systems!, Casimir forces in fluid near the critical point! and phase segregation 
in model alloys!. Many of these problems can be studied using some variation of the Ising 
model (e.g. Hopfield model of neural network, Edward- Anderson of spin glasses, etc.), or its 



magnetic ones, present often critical properties that fall in the universality class of some 
magnetic systems, the Ising model being the most simple non-trivial one. Hence, a deep 
comprehension of the general properties of the Ising model with long-range interactions is 
relevant to understand the behaviour of this kind of systems. As we will show, even the most 
simple case, i.e. the ferromagnetic model, presents non-trivial non-extensive behaviours and 
therefore it represents a good starting point to the study of more complex models. 

In this letter we consider an Ising ferrromagnet with long-range interactions, that means, 
a system described by the Hamiltonian 




Moreover, even systems not directly related with 





(1) 



with 



2 



J(r,,) = ( J > 0; a > 0) (2) 

ij 

where Vij is the distance (in crystal units) between sites i and j, and where the sum 
runs over all distinct pairs of sites on a ci-dimensional simple hypercubic lattice. The a ^ oo 
limit corresponds to the first-neighbor model. The a = limit corresponds, after a rescaling 
J ^ J/N, to the Curie- Weiss model. 

Let us introduce the sums (j)i{a) = Ji^ij)- A sufficient condition (and believed to 
be necessaryS) for the existence of the thermodynamic limit of this system is that 

0(a) = lim < (3) 

I 

Let us now take a d-dimensional hypercube of side L + 1 and N = {L + 1^, and let i = 
be the central site of the hypercube. We have that 

(j){a) = \im^(j)o{a). (4) 



Then 



L/2 L/2 ^ 



Using Euler-McLaurin sum formulalll we can approximate, for L ^ 1, 

rL/2 .L/2 1 

00(a) -J2y^ dx,...j^ dx, (,2^^2^...^^2)./. 



rL/2 

(xJ I dr r'^-^-^. 



Hence, 0o(a) shows the following asymptotic behaviour for ^ 1: 



0o(a)~ JC,(a)2" ^""/^ ^ ^ ^ (6) 



In if a = c? 



In other words: 



lim - — — a/d) (k ((^^ = CAo) for a ^ d (7) 

Ar-*oo (A^l-"/cf - 1) J2° ^ ^ ^ ^ ^ ' 
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and 



lim 



(f)oid) 



CM 



(8) 



iv^oo J2'^ In 

where Cd{a) is a continuous function of a independent of A^, with (7^(0) = 1. Therefore, 
the thermodynamic hmit is well defined for a > d (where the system presents extensive 
behaviour), while for a < d the system becomes non-extensive, the critical temperature 
becomes infinite and the standard Maxwell-Boltzmann formalism cannot be applied0lli. 
The system undergoes a second order phase transition at finite temperature for all a > d 
when0 d> 2 and for 1 < a < 2 when0 d = 1. For a ^ d^, the critical temperature shows 
the following asymptotic behaviour0: 



kBTc J (f){a) 



(9) 



We now introduce a new model that generalizes the Curie- Weiss one. Such model is 
described by the Hamiltonian: 



with 



where 



(id) 



J'in 



J{nj) 
iV*(a)2" 



(10) 
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which behaves as 



N*{a) 



1 — a/d 



N*ia) ~ <^ 



a/d-l 



ioT a/d> 1 
ioT a/d = 1 



_j_f^i-a/d < a/c/ < 1 



(12) 



(13) 



for N ^ oo. This model reduces to the Curie- Weiss one for a = and to our original model 
Eq.(|l]) (after rescaling J — > J {a/d — l)/2°') for a > d. From Eqs.(P), @) and (^ we see that 



the thermodynamic hmit of this model is well defined for all a > 0. We expect this system 
to show a phase transition at finite temperature for all a > when d >1 and for > a > 2 
when d = 1. 

The mean field theory for this model predicts a critical temperature ksT^ = J Cd{<y)2", 
which is exact for a = and for a — d~^. Hence, we conjecture the critical tempera- 
ture reproduces exactly the mean field prediction for all < a < d. This conjecture 
is difficult to verify for d > 1, since for systems with long-range interactions it is hard to 
obtain reliable numerical data for the exact critical temperature. In what follows we show 
that Ci{a) = 1 for < a < 1 and then we will test our conjecture through a Monte Carlo 
numerical simulation. 

Let us consider the d = 1 system. In this case 0o = 2 J I]n=^i Then, for a > 1 
0(a) = 2 J ({a) {C{x) is the Riemann Zeta function) and the critical temperature diverges 
a^^ ksTc ~ 2 J/ (a — 1) for Using the asymptotic behavioursEl 

n-' for Re{z) > -1 and ^ ^ 1 (14) 

n=l I — Z 

M 1 

InM (15) 

n=l ^ 

for M — i> oo, we get for a < 1 



2°4^ for < a < 1 
2 1nA^ fora = l 



(16) 



and 



for < a < 1 

Ci(«) = <i (17) 
■^::^({a) for a > 1 

For d = 1 the following necessary condition must be satisfied in order to have a finite 
critical temperature^l: 



N 



lim ^ ?7, J(?T,) = cxD. (18) 



n=l 



Using Eq.(|l^), we see that 



n J' in) ~ ; — r 

^1 ^ ' N*{a) 2^2 -a 



Hence, the critical temperature for = 1 will be finite V < a < 2. 

If we denote by m', s' and /' the energy, entropy and free energy per particle associated 
with the Hamiltonian H' i.e.. 



f'(T) = lim -^InZ', 



with Z = Tr^Si} exp{—l3 H'), 



u\T) = lim l^Tr{s.}H'expi-pH') 

iv— »oo iV 

and 

/'(T) = u'{T)-Ts'{T) 

we see that the generalized thermodynamic behaviour associated with the Hamiltonian (|1]) 
can be accommodated, for all a > 0, with the following scalings (in the limit N ^ oo and 
for T > 0): 

U{N,T) r-^ N N*u'{T*) (19) 
S{N,T) ^ N s'{T*) (20) 
F{N,T) ^ N N* f'{T*) (21) 

with T* = T/N*, as was recently conjectured by Tsallis for general systems with long- 
range interactions^. Moreover, it can be easily shown that this type of scaling preserve the 
Legendre transformation structure of the thermodynamics, even in the long-range regimeEl 
< a < d. It is also expected that the magnetization M = {J2iSi) scales as M{N,T) ~ 
Nm{T*). Therefore, the suitable plot for looking for data collapse in a numerical simulation 
will be M{N,T)/N vs T/N*. 

Let us consider the d = 1 case. We performed a Monte Carlo simulation on a chain of 
spins with Hamiltonian (0), using heat bath dynamics, for = 75, 150, 300, 600, 1200 and 
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2400. We calculated the root-mean-square of the magnetization of the system M{N, T) as a 
function of the temperature T fora = 0, 0.25, 0.5, 0.75 and 1.5. The results were averaged 
over K samples with different random number sequences ( K = 100, 50, 20, 20, 10 and 5 
for = 75, 150, 300, 600, 1200 and 2400 respectively). For every value of a we obtained 
an extrapolated magnetization curve Mq,(T), by performing a numerical extrapolation of 
M{N,T) in 1/A^ to iV ^ oo. 

In Fig.(^ we show our results for M{N,T)/N vs T/{N* 2") for a = 0.5 and 1.5. These 
curves show clearly the data collapse above mentioned. Moreover, for < a < 1 our results 
show that all curves Mq(T) fall into a single one, which coincides with the well known 
exact solution for the Curie- Weiss model {a = case), i.e. the solution of the equation 
m = tanh (m/T') (m = M/N; T' = T/{N* 2"). This last result is impressive. It does not 
only confirm our conjecture concerning the critical temperature {T^ = 1), but also shows 
that the full equation of state m = m{T) at zero magnetic field becomes independent of a 
in the non-extensive regime < a < 1, suggesting that all the thermodynamic functions are 
those predicted by the mean field theory. These results are consistent with recent Monte 
Carlo simulations of the correlation function, which reproduce the mean field behaviour in 
the same region^ < a < 1. 

In this letter we have found a new scaling for the Ising model with long range interactions 
that allows us to get a well defined thermodynamic limit for any value of a. In particular, 
for a = 0, we recover the well-known Curie- Weiss scaling, which has been vastly used in 
the context of magnetic systems. With this scaling we were able to obtain the generalized 
thermodynamic behaviour for N ^ oo (Eqs. p!9|-|21|) in the ferromagnetic case, which had 



been previously conjectured in a more general context by C. TsallisS. It is worth stressing 
that with this scaling both extensive and non-extensive behaviours can be accommodated 
in a unified and elegant formalism and the (until now almost unexplored) < a < d case 
becomes tractable. In the same way the non-extensive fully connected Ising model showed 
to be a very useful tool when suitable rescaled (Curie Weiss model), we believe that the 
model here analyzed can represent not only a useful approach but also a more realistic one 



for certain problems such as neural networks and spin glasses among others. 

On the other hand, due to the distance dependence of the interactions it becomes very 
difficult to obtain exact analytical results even in the d = 1 case. We calculated the critical 
temperature in the mean field approximation for any value of d and presented some numerical 
evidence that, (at least for d = 1) not only it reproduces the exact value in the whole non- 
extensive regime < a < c? = 1, but also the full magnetization curve M{T)/N is given by 
the mean field one. We believe that the critical temperature satisfies that property for all 
d. This conjecture is partially supported by the fact that it holds for a = and a ^ d'^. 
Moreover, since the critical exponent are those of the mean field theory both for a = 
and a ^ d'^, we conjecture that all the critical properties will reproduce the mean field 
behaviour for < a < d. Monte Carlo simulation of the correlation functionlli for d = 1 also 
support this statement. These results, although intuitive, are non-trivial and important, 
specially concerning spin glasses and biological systems (neural networks, immunology, etc) 
where a common approximation consists in to consider fully connected models instead of 
the more realistic ones with slow decaying interactions (e.g. RKKY). Our results show that 
mean field behaviour is robust again variations of the range of the interactions a whithin 
the non-extensive region, at least ioi d = 1. If our conjecture were true, this would have 
important practical implications: if you are considering systems with slow enough decaying 
interactions then you do not need sophisticated approximations, at least as far as critical 
properties are concerned! 

It would be very interesting to extend the present analysis to more general systems of 
interacting particles with long range interactions. 
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FIGURES 

FIG. 1. Monte Carlo calculations of root-mean-square of the magnetization per spin 
M{N,T)/N vs the scaled temperature T/{N* 2°^) using Hamiltonian (|l]), for a = 0.5, 1.5 and 
different values of the number of spins N in the one-dimensional lattice. In all cases the error 
bars are smaller than 0.01. The dashed lines are the N = oo extrapolation of M(N,T)/N . The 
extrapolated curves for a = 0, 0.25 and 0.75 are indistinguishable from the previous one. The solid 
line is the exact solution of the Curie- Weiss model. 
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Figure 1 - "Long-range interactions and non-extensivity.." 
S. A. Cannas and F. A. Tamarit 



